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ABSTRACT 


The  effect  of  elastic  boundary  conditions  upon  the  coupled 
bending  vibrations  and  torsional  vibration  of  a  thin  walled  beam  of  open 
cross  section  with  an  applied  axial  load  is  investigated.  Results  are 
presented  for  elastic  end  conditions  and  also  for  fixed  supports  and 
simple  supports. 
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CHAPTER  I 


INTRODUCTION 


1.1  PURPOSE  OF  THESIS 

This  thesis  reports  on  the  investigation  of  a  proposed  experi¬ 
mental  method  for  determining  the  axial  load  in  a  member  of  an  existing 
structure.  This  method  is  of  practical  importance  for  verifying  the  theo¬ 
retically  calculated  forces  in  axially  loaded  members  of  complicated 
structures  and  for  measuring  the  change  of  load  in  members  of  structures 
which  are  settling.  Plots  of  the  three  lowest  theoretical  natural  fre¬ 
quencies  versus  axial  load  are  prepared,  with  given  end  conditions  as  a 
parameter.  The  axial  load  in  a  single  span  member  may  then  be  determined 
by  comparing  its  experimentally  measured  fundamental  frequency  of  vibration 
with  the  corresponding  theoretical  plot  for  a  member  with  the  same  cross 
section  and  support  conditions. 

This  necessitates  a  detailed  investigation  into  the  effect  of 
end  conditions  upon  the  vibrations  of  a  member  of  arbitrary  cross  section 
under  various  applied  axial  loads,  which  is  the  major  topic  of  the  thesis. 

1.2  HISTORICAL  REVIEW 

The  vibration  of  a  thin  walled  open  cross  sectioned  member  under 
an  eccentrically  applied  axial  load  were  discussed  by  V.Z.  Vlasov  (1,2)  in 
1940.  He  derived  the  differential  equations  for  coupled  vibrations  and 
solved  them  for  the  case  of  simple  supports  which  in  this  thesis  means  that 
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the  member  is  not  allowed  to  deflect  perpendicular  to  the  longitudinal 
axis  or  to  rotate  about  the  longitudinal  axis.  Vlasov's  solution  em¬ 
ployed  the  method  of  separation  of  variables  and  assumed  the  deflected 
form  to  be  a  sine  function.  Timoshenko  (3)  considered  the  problem  of 
finding  the  natural  frequencies  for  the  coupled  torsional  and  bending 
vibrations  for  a  channel  cross  section  having  simple  supports.  In  1940 
Garland  (4)  used  the  method  of  Rayleigh-Ritz  to  determine  the  frequencies 
and  amplitude  ratios  of  a  cantilever  beam  with  channel  cross  section.  He 
verified  his  results  by  experiment.  Frederhoffer  (5)  solved  the  general 
case  for  a  symmetrical  cross  section.  A  solution  for  an  arbitrary  cross 
section  using  energy  methods  was  presented  by  Karyakin  (6). 

In  1954  Gere  (7)  analysed  the  free  torsional  vibrations  of  bars 
of  thin  walled  open  cross  sections  for  which  the  shear  center  and  the 
centroid  coincide.  He  investigated  the  effect  of  warping  on  the  frequency 
of  single  span  members  with  various  end  conditions.  Gere  and  Lin  (8) 
solved  the  differential  equations  for  coupled  free  vibrations  of  beams  of 
non  symmetric  open  cross  section.  They  considered  three  different  end 
conditions:  simple  supports,  fixed  ends  and  a  cantilever.  They  solved 
their  derived  equations  by  an  exact  method,  using  a  computer  to  do  the  calcu¬ 
lations.  They  also  used  the  Rayleigh-Ritz  method  to  obtain  an  approximate 
solution. 

Aggarwall  and  Cranch  (9)  further  extended  the  theory  of  coupled 
bending  torsional  vibrations  by  including  the  effects  of  rotatory  inertia 
and  transverse  shear  deformation.  Their  results  agree  with  the  Timoshenko 
theory  for  the  low  frequencies  and  give  more  satisfactory  results  for  high 
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frequencies  since  it  predicts  finite  velocities  of  wave  propogation  for 
high  frequencies  where  as  the  Timoshenko  theory  predicts  infinite  velocities. 

Didrikson  (10)  applied  the  computer  method  of  Gere  and  Lin  to 
the  problem  of  thin  walled  members  of  open  cross  section  with  an  axial 
load  applied  at  the  centroid  of  the  cross  section.  Didrikson's  experi¬ 
mentally  determined  frequencies  were  less  than  those  calculated  assuming 
fixed  ends  but  greater  than  those  calculated  assuming  simple  supports. 

This  indicated  that  the  actual  support  conditions  were  between  the  fixed 
and  simple  cases,  and  that  the  theory  should  be  modified  to  account  for 
elastic  end  conditions.  Rogers  (11)  discussed  elastic  end  conditions  for 
bending  vibrations. 

1.3  OUTLINE  OF  INVESTIGATION 

The  vibrations  considered  are  torsional  as  well  as  flexural  since 
the  member  was  of  unsymmetric  thin  walled  open  cross  section.  Torsion  of 
a  beam  generally  produces  warping  of  the  cross  section.  Warping  is  the 
movement  in  the  longitudinal  direction  of  a  point  out  of  the  plane  perpendicu¬ 
lar  to  the  longitudinal  axis  of  the  member.  If  a  member  has  torque  applied 
only  at  the  ends  and  all  cross  sections  are  free  to  warp  then  the  member  is 
said  to  be  in  pure  torsion.  A  member  is  under  nonuniform  torsion  if  the 
torque  varies  along  the  member  or  if  some  of  the  cross  sections  are  restrained 
from  warping.  In  nonuniform  torsion  the  warping  varies  along  the  length 
of  the  member  so  that  axial  stress  accompanies  the  torsional  shear  stress. 

The  torsional  and  flexural  vibrations  are  coupled  when  the  axis 
joining  the  shear  center  of  the  cross  sections  of  the  member  does  not  coin- 
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cide  with  the  longitudinal  axis.  The  axis  through  the  shear  center  is 
the  only  axis  of  the  member  where  an  applied  transverse  load  will  produce 
flexural  deformations  only.  Since  the  flexural  and  torsional  frequencies 
are  influenced  by  the  end  supports  of  the  member,  support  conditions  for 
a  vibrating  member  are  discussed  in  detail  in  the  following  chapter. 


CHAPTER  II 


THEORY 

This  chapter  outlines  the  derivation  of  the  differential 
equations  describing  the  free  coupled  bending  torsional  vibrations  of  a 
thin  walled,  open  cross  section  member  under  a  concentric  axial  load. 

These  differential  equations  are  solved  for  the  lower  natural  modes  of 
vibration  employing  the  method  used  by  Gere  and  Lin  (8)  when  they  solved 
the  problem  for  zero  axial  load.  The  boundary  conditions  considered  in 
this  solution  are  elastic.  Ordinary  differential  equations  are  obtained 
from  the  partial  differential  equations  by  the  method  of  separation  of 
variables.  Rotatory  inertia  and  shear  deformation  are  neglected  since 
only  the  lower  fundamental  frequencies  are  investigated. 

In  the  derivation  of  the  differential  equations  it  is  assumed 
that  the  member  has  thin  walls,  is  long  and  slender  with  a  constant  cross 
section,  and  is  composed  of  a  homogeneous  material.  It  is  further  assumed 
that  the  angle  of  rotation  and  the  transverse  deflections  are  small. 
Bending  is  taken  about  the  centroidal  axes  and  twisting  about  the  longi- 

c 

tudinal  axis  through  the  shear  center.  The  thrust  is  assumed  to  act  along 
a  line  parallel  to  the  longitudinal  axis  of  the  member. 

2.1  DERIVATION  OF  GOVERNING  EQUATIONS 

The  coordinate  system  used  in  the  derivation  of  the  equations 
is  shown  in  figures  2.1  and  2.2.  The  y  and  z  axes  have  their  origin  at 
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FIGURE  2.1  COORDINATE  SYSTEM 


Figure  2.2  ORIENTATION  OF  PRINCIPAL  AXES 
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the  shear  center  and  are  parallel  to  the  n  and  c  axes  respectively.  The 
n  and  c  axes  are  the  principal  axes  of  the  cross  section  with  their  origin 
at  the  center  of  gravity.  The  x  axis  is  the  longitudinal  axis  of  the 
member  through  the  shear  center  of  the  section.  The  deflection  in  the  y 
and  z  directions  are  v  and  w  respectively  with  the  angle  of  twist  about 
the  longitudinal  axis  being  denoted  as  <J>.  The  distances  from  the  shear 
center  to  the  center  of  gravity  are  c^  and  cz  in  the  y  and  z  directions 
respectively. 

The  bending  of  the  member  about  the  z  and  y  axes  by  a  trans¬ 
verse  load  is  expressed  by 


El 


C 


El  =  -  M 

n  d7  y 


2.1 


2.2 


The  theory  of  nonuniform  torsion  of  a  member  with  a  thin  walled 
open  cross  section  is  here  discussed  as  it  is  necessary  for  derivation  of 
the  expression  for  torque  on  a  cross  section. 

When  a  bar  is  subjected  to  pure  torsion  the  torque  is  given 

by 


M 


t 


Gc0 


2.3 


where 


G  is  the  torsional  rigidity  of  the  bar 


6  is  the  constant  angle  of  twist  per  unit  length 
c  is  the  torsion  constant 
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The  value  of  the  torsion  constant  for  a  cross  section  which  consists  of  n 
thin  bars  of  different  thickness  is  given  by 


where  is  the  length  of  the  middle  line  of  the  ith  member 

t.  is  the  thickness  of  the  ith  member 


The  warping  displacement  u  for  pure  torsion  is  given  by 
Timoshenko  (3)  as 


u  =  e(w„  -  w) 
s  s 

The  quantity  w$  is  called  the  warping  function  and  is  defined  by 

s 

w  =  /  r  ds 
s  Q  z 

where  rt  is  the  perpendicular  distance  from  the  tangent  at  a  point  M 
on  the  midline  of  the  member  cross  section  to  the  axis  of 
rotation  as  shown  in  figure  2.2, 

s  is  the  distance  along  the  middle  of  the  member  to  the  point  M 
under  consideration. 

The  quantity  wsis  the  average  value  of  w  and  is  therefore  defined  as 
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ws  =  F  l  ws  ds 


where  b  is  the  total  length  of  the  midline  of  the  cross  section. 

For  nonuniform  torsion  the  constant  angle  of  twist,  0,  is  re¬ 
placed  by  the  rate  of  change  of  rotation  about  the  longitudinal  axis, 
^  ,  in  the  formula  for  warping  displacement. 


u  ■  af  (ws  -  V 


2.4 


The  strain  due  to  the  warping  displacement  u  is 


£x  =  dx  =  ^2  (ws  '  ws> 


2.5 


Assuming  the  material  obeys  Hookes  Law  the  longitudinal  stress  due  to 
warping  is 


r  dl  /-  \ 

=  E  —f  (w  -  w  ) 
dx^  s  s 


2.6 


This  longitudinal  stress  produces  a  shearing  stress  x  which  produces  a 
torque  in  the  cross  section.  The  shearing  stress  is  given  by  Timoshenko  (3) 
to  be 


T 


E  d34> 
1  dx3 


s  _ 

/  (w  -  w  )  t  ds 
0  s  5 


Thus  the  torque  due  to  warping  of  the  cross  section  is 


where 


cw  is  the  warping  function  of  the  cross  section,  and  is  given 

by 


* 


c  = 
w 


b 

J  K 
o 


ws) 


t  ds 


The  total  torque  acting  on  a  bar  subjected  to  nonuniform  torsion 
is  obtained  by  adding  equation  2.3  to  equation  2.7 


=  Gc  ^  -  Ec  ^4 
dx  w  dx3 


2.8 


Differentiating  equation  2.8  with  respect  to  x  gives  the  equation  of  non- 
uniform  torsion  under  static  torque 


-  W 


t 


Gc  h  -  Ec  A 
dx2  w  dx4 
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A  member  of  thin  walled  open  cross  section  under  arbitrary  lateral 
loading  deforms  according  to  equations  2.1,  2.2,  and  2.8. 

The  effect  of  an  axial  load  upon  the  bending  and  torsion  equation 
is  now  shown.  Consider  the  point  M  on  the  mid-surface  of  the  member  cross 
section  as  shown  in  figure  2.3.  After  the  axial  load  is  applied  the  member 
deforms  with  the  point  M  moving  vs  and  w$  in  the  y  and  z  directions  respectively 
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Figure  2.3  MEMBER  OF  ARBITRARY  UNSYMMETRICAL  OPEN 
CROSS  SECTION  SUBJECTED  TO  AXIAL  THRUST 
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vs  =  V  -  <f)Z 
ws  =  w  +  4>y 

The  axial  load  on  the  member  is  positive  in  compression  and 
is  given  by 


b 

P  =  /  at  ds 
0 


where  t  is  the  thickness. 

Let  pz  denote  the  difference  of  the  projection  of  at  ds  on  the  z  axis  and 
p  the  difference  of  the  projection  upon  the  y  axis.  Thus 


where 


p  dx  ds  =  at  ds  — 
z  pz 


p  ds  dx  =  at  ds  — 
y  Py 


2.10 


p  and  p  are  the  radii  of  curvature  of  the  proejction  of  v  (x) 

z  y  s 

and  w  (x) . 
s 


For  small  deflections  the  radii  of  curvature  are  defined  by 


,2 

1  d  ws 


1 


P. 


dx2  p\ 


d2v. 


dx' 


Substituting  for  the  curvature  and  dividing  by  dx  ds  gives 


2  2 
.  rd  w  .  dd)-| 

P  =  -  at  [— T  +  y  —f] 

1  dxz  dxz 


teb 
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p  =  -  at 


These  loads  give  rise  to  a  torque  dW^  about  the  x  axis.  The 
value  of  the  torque  is  the  load  times  its  moment  arm  from  the  x  axis 
and  pz  times  its  moment  arm.  Therefore 


dW,  =  p  v  ds  -  p  z  ds 
t  vzJ  Ky 

dW.  =  -  at  [^-4  +  y  ^-4]  yds 
L  dx^  dx^ 

+  at  [^-4  -  z  ^-4]  zds 
dx  dx^ 


Integrating  dW^  over  the  cross  section  and  noting  that 


/  ytds  =  Ac  /  ztds  =  Ac  /  zytds  =  0 


A 


A 


A 


/  y^tds  =  I  /  z^tds  =  I, 


A 


A 


the  total  torque  (W^)  due  to  the  axial  load  is  obtained 


(W. )  =  a  ^-4  c  A  -  a  I  -^-4  -  a  ^-4  c  A  -  a  I  — ? 

4  p  dx2  z  n  dx2  dx2  y  ?  ^2 


Q  § 

dx‘ 


2  2  PI  2 

(W  )  =  ^  pc  .  dfw  p  .  JM 

(Vp  dx2  Pcz  Hv2  Pcy  A  ,  2 


dx 


dx 


2.11 


The  bending  moments  caused  by  the  axial  load  p  are 
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(Mc)p  =  P(v  +  cz  d>)  2.12 

(Mn)p  =  P(w  -  cy  <D)  2.13 


The  bending  equations  for  a  member  with  an  arbitrary  lateral 
load  and  an  axial  thrust  are  obtained  by  combining  equation  2.12  with  2.1 
and  2.13  with  2.2.  The  resulting  equations  are 

j2 

El  =  _  m  p(v  _  ^ 

g  dx  ^  z 


,2 

EI  =  -  M  -  P(w  +  c  <(>) 

n  dx2  n  y 


These  equations  are  expressed  in  terms  of  the  loads  in  the  y  and  z  directions 
by  differentiating  twice  with  respect  to  x  to  obtain 


4  dx^  y  z  - ~c 


2.14 


dx 


dx 


EI  ^4  =  W  -  P  £$.  c 
n  dx4  z  h*2  y  -.2 


2.15 


dx 


dx 


The  torque  equation  obtained  by  combining  equation  2.11  with  2.9  is 
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-W 


t 


Ec  4  + 
w  dx4 


(Gc  - 


d2w 

dx2 


+  Pc. 


d2v 


2.16 


For  free  vibrations  the  only  loads  acting  on  the  member  other 
than  the  axial  load  are  the  loads  due  to  inertia  forces.  Therefore  the 
components  of  the  loading  in  the  transverse  direction  are 

W  =  -  mA  (v  -  c  <j>) 
y  8t  z 

a2 

W  =  -  mA  — p  (w  +  c  <J>) 

at  y 


The  inertial  torque  is 


2  2  2 

(W.  )T  =  -  I  +  mA  (v  -  4>c  )c  -  mA  — 

L  1  p  at  at  at 


2  (W  +  Cy^  Cy 


a2v 


a2w> 


=  .mI  ^+mA(c  z^-c 

at  z  at  y  at 


When  the  inertial  loads  are  substituted  into  equations  2.14, 
2.15  and  2.16  the  differential  equations  describing  the  free  vibrations 
of  a  thin  walled  member  of  open  cross  section  under  an  applied  axial  load 
are  obtained.  They  are 

EI  A  +  P  4  +  "A  4  .  Pc  4  .  mAc  ifi  .  0 

^  ax  ax  at  z  axz  z  at 


2.17 
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El  2JJ.  +  p  +  mA  sfw  +  p  A  +  mAc 
n  ^4  ^2  ^2  y  av2  y 


3x 


9x 


at 


OL  =  0 
2  u 
at 


2.18 


4  PI  2  2  2  2 

r„  a  <j)  ,  r  o  3i  ,  _T  ai  n_  a  v  „n_  a  v 

^cw  4  +  L  /\  "  Gc]  2  +  mIo  2  ~  ^cz  2  "  m^Cz  2 

w  ax4  A  ax^  u  at*  z  ax^  2  at* 


+  Pc  +  mAc  =  0 

y  9x^  y 


at 


2.19 


The  above  equations  do  not  consider  the  effects  of  rotatory 
inertia  or  shear  deflection.  Both  effects  were  considered  by  Aggarwal 
and  Cranch  (9)  for  a  member  with  zero  axial  load.  They  stated  that  ro¬ 
tatory  inertia  and  shear  deflections  have  a  large  effect  only  upon  the 
high  vibrational  frequencies.  In  this  thesis  only  the  three  lowest  natu¬ 
ral  frequencies  are  considered  so  the  terms  accounting  for  rotatory  inertia 
and  shear  deflections  are  neglected.  The  equations  derived  are  similar  to 
those  given  by  Vlasov  (2). 

2.2  BOUNDARY  CONDITIONS 

Before  considering  the  case  of  elastic  boundary  conditions,  the 
limiting  cases  of  simple  supports,  fixed  supports  and  free  ends  will  be 
di scussed . 

Simple  supports  do  not  allow  any  deflection  of  the  member  at  the 
supports,  but  they  do  not  offer  any  resistance  to  rotation  of  the  member 
about  an  axis  perpendicular  to  the  longitudinal  axis.  A  simple  support 
provides  restraint  against  rotation  about  the  longitudinal  axis.  Warping 
of  the  cross  section  is  not  restrained,  which  makes  the  stress  in  the  longi¬ 
tudinal  direction  due  to  warping  equal  to  zero.  Considering  the  above  de- 


finition  of  simple  supports  and  equation  2.6  for  the  stress  due  to  re¬ 
strained  warping,  the  boundary  conditions  for  simple  supports  are  seen 
to  be 
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v  =  w  =  =  0 


2  2  2 

3  v  _  3  w  _  a  <£> 

2  ~  2  "  2 
3xd  3yd  9x 


0 


A  fixed  support  is  one  which  clamps  the  member  rigidly.  There 
is  complete  resistance  to  translation  of  the  member  and  to  any  rotation 
about  an  axis  perpendicular  to  the  longitudinal  axis  of  the  member.  A 
fixed  support  does  not  allow  any  rotation  about  the  longitudinal  axis  nor 
any  warping  of  the  cross  section.  Equation  2.4  is  zero  when  the  warping 
displacements  are  zero.  Thus  the  boundary  conditions  for  a  fixed  support 
are 


v  =  w  =  (J>  =  0 

3v  _  3w  _  3<j>  _  n 
3x  ~  3x  9x 

When  the  support  condition  at  one  end  of  the  member  is  free 
there  is  no  resistance  to  translation  or  rotation.  This  means  that  the 
shear  at  the  end  and  the  bending  moment  at  the  end  must  be  zero.  Also, 
the  longitudinal  stress  due  to  warping  and  the  total  torque  must  be  zero. 
From  equations  2.6  and  2.8,  and  the  above  requirements  the  boundary  condi- 


18 


tions  for  a  free  end  are 


2  2  2 

9  v  _  9  w  3'  4)  _  A 

2  2  2  U 
9x  9x  9x 


3 


=  0 


Gc 


3(j) 

9x 


A  member  on  elastic  supports  may  be  considered  to  be  on  a  set 
of  supporting  springs.  One  spring  resists  transverse  deflections  of  the 
member  and  another  resists  rotation.  The  force  in  the  spring  resisting 
deflections  must  be  equal  to  the  shearing  force  at  the  end  of  the  member, 
and  the  moment  supplied  by  the  spring  resisting  rotation  of  the  end  must 
equal  the  bending  moment  at  that  end.  A  third  spring  resists  the  displace¬ 
ments  due  to  warping  and  a  fourth  resists  rotation  about  the  longitudinal 
axis  of  the  member.  Therefore,  the  force  applied  by  the  third  spring 
equals  the  resultant  of  the  longitudinal  stress  in  the  beam  due  to  warping 
of  the  cross  section,  and  the  torque  in  the  fourth  spring  equals  the 
torque  at  the  end  of  the  beam.  These  conditions  are  expressed  mathematically 
as 


V  =  k-|U 


m  =  a-,0 
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where 


u  and  0 
member. 


y  (ws  -  ws>  =  k2  If  (ws  -  ws) 


G  c  +  Ec...  =  3o4> 


8x 


'W3X3  M2 


u  is  the  vectoral  displacement  of  the  member 
9  is  the  angle  of  rotation  of  the  member  about  an  axis 
perpendicular  to  the  longitudinal  axis 
k-| ,  3-| »  ancl  32  are  sPr‘'n9  constants 

can  be  broken  into  components  parallel  to  the  given  axes  of  the 
When  this  is  done  the  boundary  conditions  become 


33v 

El  +  k,v  =  0 

c  3xJ  1 


El  ^4+  k,w  =  0 
n  3x3  1 


EI  +  B,  Iv  =  0 
n  3x2  1  3x 


+  k  M.  =  o 
3x2+k23x 


Gc  +  Ec  2-f  -  B,  <t>  =  0 
3x  w  3x3  2 


These  equations  incorporate  the  three  boundary  conditions  pre- 
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viously  discussed  as  special  cases.  For  the  case  of  simple  supports  the 
values  of  the  spring  constants  k-j ,  3-j »  k^,  and  82  are  infinity,  zero,  zero 
and  infinity  respectively.  For  fixed  ends  the  values  of  the  supporting 
spring  constants  are  all  equal  to  infinity,  and  for  free  ends  the  spring 
constant  values  are  all  zero. 

2.3  SEPARATION  OF  VARIABLES 


The  governing  partial  differential  equations  2.17,  2.18,  and  2.19 


can  be  reduced  to  a  set  of  ordinary  differential  equations  by  the  appli¬ 
cation  of  the  method  of  separation  of  variables.  Let 


v ( x , t )  =  V(x)  T (t) 


w(x,t)  =  W(x)  T(t) 


<j)(x,t)  =  $(x)  T(t) 


When  these  expressions  are  substituted  into  the  equations  they  yield  the 
following  ordinary  differential  equations 


2.20 


2.21 


21 


2.23 


2 

where  pn  is  the  constant  of  separation. 

The  solution  of  the  time  equation  2.20  is 


T  =  A  cos  pt+Bsinpt 

m  m 


2 

This  shows  that  the  constant  of  separation  pn  is  the  square  of  the  circu¬ 


lar  frequency  of  vibration. 

2.4  SOLUTION  OF  THE  EQUATIONS  WITH  ELASTIC  BOUNDARY  CONDITIONS 

The  method  of  solution  of  the  equations  will  be  outlined  for 
elastic  boundary  conditions.  This  method  can  be  applied  to  the  three 
limiting  cases  by  the  proper  choice  of  the  supporting  spring  constants. 

The  ordinary  differential  equations  2.21,  2.22,  and  2.23  can  be 
rewritten  using  the  operator  D  =  ^  in  the  form 


(El  D4  +  PD2  -  mApn2)V  +  (Pcz  D2  +  mApp2  cz)<D  =  0 


2.24 


( EI^  D4  +  PD2  -  mApn2)W  +  (Pcy  D2  -  mApn2  c^)$  =  0  2.25 


+  [Pc  D2  -  mAp  2  c  ]W  =  0 
L  y  Kn  yJ 


2.26 


90  r 

■ 

V 
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W  and  0  are  eliminated  from  these  equations  to  obtain  the  following 
twelfth  order  differential  equation 


E3c  II  19  c2  c  I  II  Gel  I  ln 
— y-?--  --  D  2V  +  -I-S-  {P— y— +  -^] - y^-}  D  °V 

mVl0  mV  l0  A  l0 


2  2 

p  t  cl  II  CD  c  I 

{_  ID _  rJU.  +  JLJli  +  _EP_ 

m2A2  L  !0  A  m3A2  ‘0 

p„2e 


[p [JL  +  _£  .  J_  (c  2i  +  c  2I  )] 
inl„  A  lQ  ^cz  cy  VJ 


Gel  0  p  E  c  I  ,  0  0 

-r-^]}  D8V  +  {-  -2—  [2P[-p  +  -P~-  ±  (c  2I  +  c  2I  )] 
‘o  m2A  *0  fl  J0  z  n  y  C 


Gel  d2  PI  r  /-  p„^p  3PI  9r 

-y6- ]  +  4?  -  r]}D  v  +  {-  -f~ 

X0  nfV  x0  x0  ni  A  A0A  X0 


+  ^^+/-17(cz21  +cy21  )]}°4V 


0 


0 


^n  r  P  Gc-i  n2w  „  6  p  _  n 

+  — [l/-y][)''-Pn  T^-° 


2.27 


The  system  of  equations  can  be  used  to  produce  the  same  equation 
with  W  or  $  interchanged  for  V.  Solutions  to  the  twelfth  order  differential 
equations  can  be  expressed  in  the  following  forms 


V 


A  e 
m 


V 


V 


3> 


12 


l 

m=l 


C  e 


m 


r  x 
m 


The  characteristic  equation  of  the  problem  is  derived  by  substi- 


X  ' 
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tuting  the  solution  for  V  into  equation  2.27.  It  is 


E3c  II 


3  A  2  T 
m  A  I 


0 


c  I  II  Gel  I  ln 

_  w  p  .  C  Tin  C  n-j  JO 

32LrLI  AJ"T  'rr 

mV  l0  A  *0 


E  {P[41tX3 


2  2 

p„  E^  cl  II 


+  {.  in _ r_w_£.  +  + 

1  2.2  L  In  A  J 


EP 


_  t  „  i7[p[r+i-r(cA  +  cv\)] 

m  A  A0  "  mV  l0  A  x0  z  n  y  c 


Gel  o  d2  PI  r  p  3E  c  I  ,  0  0 

-y^]}r8  +  f-JT  ^  -  f]  -  -2“  [2PCf  +  /  '  f  (cz  'n  +  C,  V] 

A0  mV  Ai0  X0  rri  A  X0  A  *0  Z  n  y  C 


GcI-  6  .  rpn  E  -c-  1 


iV+ 


p  3PI 

in_  r _ £ 

2  LI  A 
i/ A  10A 


2Gc]}r4  + 


I 


0 


m  L  I 


3PI 
[— * 


0 


f]r2 

*0 


P6^ 

"  !o 


=  0  2.28 


? 

If  the  square  of  the  frequency  pn  is  known  then  the  equation 
can  be  solved  for  the  roots  which  will  be  in  positive  and  negative  pairs 
since  only  even  powers  of  r  appear.  Therefore,  the  solution  may  be  written 
as 

r-,x  -r-,x  r?x  -r?x  -rfix 

V  =  A-je  +  A^e  +  A3e  +  A4e  A-j  2e 


r,x  -r-jX  -rfix 

W  =  B-je  +  B^e  +  ........................  +  B^e 


r-jX  -r-jX  -rgx 

^  =  C  i  e  +  C2e  +  ........................  +  C-j  2e 


The  arbitrary  constants  B.  and  can  be  written  in  terms  of  the 


2*19 
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constants  A.  by  substituting  the  solution  for  V,  W,  and  $  into  equations 
2.24  and  2.25  to  obtain  the  following  equations 


[El  r' 


+  Pr-j  -  mAp 


2  rlx 
^]A-|  e  1 


-  [Pczr] 


\ipn 


'  Cz ] c i e 


rlx 


42  9  -r I x  2  2  ~rlx 

[EI^rl  +  Prl  "  mAPn  ]A2e  "  [pc2r]  *  mApn  cz]C2e 


+  [EIrr64  +  Pr62  "  mApn^A12e  6  '  [pczr]2  *  mAPn2cz3c12e  6=0 


4  2  2  rlX  2  2  rlx 

[El^r 4  +  P-j r]  -  mApn  +  [Pc  ^  -  mApn  c  ]C.,e 


+  [EInri4  +  piri2  '  mApn2]B2e  1  +  [Pc  -  mApp2c  ]C2e  1 


-r,x 


o  •  o  •  e 


+  [EInr64  +  P^2  -  tnApn2]B12e  6  +  [Pc  r&2  -  mAp/c  ]C]2e  6  =  0 


rl x  -r,x  r„x 
Since  e  ,  e  ,  e 


_r6x 

e  are  linearly  independent 


the  coefficients  of  each  must  be  zero.  This  enables  the  constants  and 


B.j  to  be  expressed  in  terms  of  A...  Let 


E  Irrl  P  2  2 

r  r  ^  1 _  +  '  .  r  c  -  n  c 

,  L1  u2  mA  mA  1  P' 

,"i  =  a7  =  r2  =  - 


n 


rl2  -  pn2) 
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r  r  EIC  4  P  2  2 

,  _  l3  _  l4  _  mA  r2  mA  r2  “  p* 

^2  '  A  '  A  -  (_L  r  2  2} 

J  4  VmA  r2  pn  ‘ 


n 


B-,  B, 


X1  ~  C- 


El 


^  -  Pn') 


n  4  ,  p  2 


ri  +  -  P. 


mA  '1  mA  ' 1 


B„  B 


x2  C0  C 


El 


Pr  c 

r  r  2  2 

cy(~sr  -  pn  > 


1  „  4 


r0  + 


mA  '2  '  mA  r2  *  p 


El. 


Put  y-j  =  X]  = 


h. 

A1  C1 


^2  ^2 

a2  c2 


4  .  Prl 


mA  '  1 


rn  +  — r« - p 


mA 


EIn  4  PV 
mA  rl  mA  pr 


C3  B3 


C4  B4 


Y2  ~  1 2X2  ~  c  c 


EIC  4  Pr2 
mA  r2  mA 


-  P 


n 


El  A  Pr0  0 

n  r„4  +  -4-  -  P  2 


mA 


mA 


Using  the  ratios  y.  and  ib.  the  coefficients  B.  and  C.  can  be 

3  'll  li 


written  in  terms  of  the  coefficients  A^  since 


B .  =  y.A. 

i  'ii 


Ci  =  *iAi 


The  values  of  the  roots  r  will  not  all  be  positive;  therefore. 
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some  of  the  exponents  will  be  imaginary  and  the  exponentials  may  be 
combined  to  form  sines  and  cosines.  The  solution  will  take  the  following 
form,  assuming  the  second  root  is  negative 


r  x  "rlx  1  1 

,,  A  Q  1  +  A0e  +  A,  cosr0x  +  A.  sinr0x  + 

v  =  Aie  2  3  2  4  2 


rlX  "rl^  1  1 

W  =  y-| A-j e  +  y-]A2e  +  y2A3  cosr2x  +  ^2^4  s‘>nr2x  +  .  2.29 


r-jX 


"rlX  ,  ,  n  1 


1  ,, 


$  =  A-j e  +  ip-j A2e  +  ij^Ag  cosr2x  +  ij^A^  sinr2x  + 


where  =  A3  +  A^ 


1 


a4‘  =  /-  1  (A3  -  a4) 


The  unknowns  of  the  problem  as  stated  are  pn  with  a  corresponding 
twelve  values  of  r,  and  for  each  pair  pn  and  r  these  are  twelve  constants 
A..  There  are  thirteen  equations  derived  to  solve  for  the  unknowns.  These 
equations  are  the  boundary  equations  at  each  end  of  the  member  and  the 
characteristic  equation  2.28.  If  a  complete  solution  to  the  problem  were 

v  V'\ 

desired  then  it  could  be  found  in  terms  of  one  of  the  constants  A^ ,  which 
in  turn  could  be  found  by  the  use  of  a  given  set  of  initial  conditions  and 
equation  2.20.  Of  main  interest  in  this  thesis  is  a  solution  for  the 
natural  frequencies  of  the  equation  which  do  not  depend  upon  the  initial 
conditions  of  the  problem,  but  do  depend  upon  the  boundary  conditions.  The 
following  symmetric  elastic  boundary  conditions  are  considered 


' 
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EIc  77  *  kiv  ■  0 
s  dx 

at  x  =  0  and  x  =  £ 

3 

d  W 

El  ^-4  ♦  k,W  =  0 
n  dx3  1 

at  x  =  0  and  x  =  £ 

El  dV  -  n 

EI?  dx2  61  dx  0 

at  x  =  0 

El  d^V.  +  o  dV  -  n 
e  dx2  61  dx  -  0 

at  x  =  £ 

El  dfw  .  dW  .  Q 

EIn  dx2  dx  -  0 

at  x  =  0 

El  ^  +  6,  7p  =  0 
n  dx2  1  dx 

at  x  =  £ 

d2i  .  dt  .  n 
dx2  2dx 

at  x  =  0  and  x  =  £ 

Gc  37  +  Ecw  -  62*  ■  0 

at  x  =  0 

Gc  37  +  Ecw  +  62*  ■  0 

dx 

at  x  =  £ 

Twelve  homogeneous  equations  are  obtained  when  the  solution 
for  V,  W,  and  $  are  substituted  into  the  above  boundary  conditions.  For 
a  nontrivial  solution  the  determinant  of  the  coefficients  must  be  zero. 
This  determinant  will  be  called  the  boundary  condition  determinant.  When 
the  solutions  to  V,  W  and  $  are  in  the  form  of  equations  2.29  the  boundary 
condition  determinant  will  be  of  the  following  form 


A 


■•iu 


o 

II 
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If  a  natural  frequency  pn  is  known  the  characteristic  equation 
can  be  solved  for  the  roots  r. .  When  these  roots  are  substituted  into 
the  boundary  condition  determinant  its  value  will  be  equal  to  zero.  The 
use  of  the  boundary  condition  determinant  and  the  characteristic  equation 
suggests  a  trial  and  error  method  of  solving  for  the  natural  frequencies 
of  a  member.  First  a  value  of  thrust  is  selected  and  a  value  of  frequency 
is  assumed.  These  are  substituted  into  the  characteristic  equation,  which 
is  then  solved  for  its  roots.  These  roots  are  substituted  into  the  boundary 
condition  determinant,  and  the  determinant  is  evaluated.  If  the  value  is 
not  equal  to  zero  then  another  frequency  is  assumed  until  one  is  found 
which  makes  the  determinant  equal  to  zero.  This  frequency  is  a  natural 
frequency  of  the  member.  A  digital  computer  was  employed  to  facilitate 
the  calculations. 

This  method  of  solution  was  used  to  solve  the  specific  problem 
of  an  angle  member  with  an  unequal  leg  cross  section  with  the  bending 
moments  at  the  ends  resisted  elastically  and  the  rest  of  the  supports  the 
same  as  fixed  ends.  The  boundary  conditions  were  taken  as  symmetrical. 

2.5  SOLUTION  FOR  MEMBER  OF  UNEQUAL  LEG  ANGLE 

For  an  angle  member  with  an  unequal  leg  cross  section  the  value 
of  the  warping  constant  cw  is  zero,  since  the  cross  section  is  of  thin 
rectangular  elements  with  a  common  point  of  intersection  as  stated  by 
Timoshenko  and  Gere  (12).  When  the  warping  constant  is  zero  the  governing 
differential  equation  is  reduced  to  tenth  order,  so  the  characteristic 
equation  is  reduced  to  tenth  degree.  The  characteristic  equation  is 
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2.30 


The  values  of  ^  and  y  are  not  affected  by  the  value  of  c  . 

The  size  of  the  boundary  condition  determinant  is  reduced  by  the  fact  that 
cw  is  zero.  Since  the  characteristic  equation  is  of  tenth  degree  there 
are  only  ten  roots;  therefore,  two  columns  are  dropped  from  the  determinant. 
When  c  ,  is  zero  it  means  that  (w  -  w  )  is  equal  to  zero;  therefore,  the 
equations  describing  the  warping  displacements  of  the  cross  section  are 
identically  equal  to  zero  and  cannot  be  used  as  boundary  conditions.  This 
reduces  the  number  of  rows  in  the  boundary  condition  determinant  by  two. 

When  cw  is  equal  to  zero  as  it  is  with  an  unequal  leg  angle  cross  section 
the  boundary  condition  determinant  is  reduced  to  ten  rows  and  ten  columns. 

The  boundary  conditions  used  were  those  thought  to  most  closely 
approximate  the  boundary  conditions  of  the  experimental  work  done  in  this 
thesis.  Since  there  was  no  translation  at  the  supports  the  value  of  k-j 
was  taken  to  be  infinity.  The  supports  were  assumed  to  allow  no  rotation 
about  the  longitudinal  axis  so  the  value  of  3 2  was  set  equal  to  infinity. 
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Since  comparison  of  experimental  results  with  theoretical  results  for 
fixed  supports  showed  that  the  supports  were  elastic,  the  value  of  3-j 
was  taken  non  zero  and  finite.  The  ten  resulting  boundary  conditions 
used  were 


V  =  W  =  $  =  0 


at  x  =  0  and  x  =  i 
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The  determinant  which  results  from  using  these  boundary  conditions  is 
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The  trial  and  error  solution  outlined  in  the  previous  section 
was  used  to  solve  2.30  and  2.32  for  the  three  lowest  natural  frequencies. 
An  IBM  360  computer  was  programmed  in  the  Fortran  IV  language  to  carry 
out  the  necessary  computations.  The  data  read  into  the  program  were 
the  values  of  thrust  and  frequency  squared.  The  program  was  run  once 
for  each  value  of  3.  The  results  were  plotted  as  thrust  versus  the 
square  of  the  natural  frequency,  with  the  value  of  6  used  as  a  parameter. 
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CHAPTER  III 


EXPERIMENTAL  APPARATUS  AND  PROCEDURE 

The  experiments  conducted  were  designed  to  test  the  effect  of 
different  support  conditions  upon  the  natural  frequencies  of  a  single  span 
thin  walled  member  of  open  cross  section  with  an  applied  axial  load.  The 
value  of  the  axial  load  was  varied  from  near  the  compressive  buckling  load 
of  the  member  to  ten  thousand  pounds  tension,  with  the  three  lowest  natural 
frequencies  being  determined  for  each  load.  The  test  member  was  an  un¬ 
equal  leg  angle  which  was  bolted  by  various  methods  to  a  fixed  support  at 
one  end  and  to  a  sliding  support  at  the  other.  A  hydraulic  pump  and  ram 
system  was  used  to  move  the  sliding  support  and  thus  applied  the  axial 
load  to  the  member.  The  strain  caused  by  the  axial  load  was  measured 
by  electrical  resistance  strain  gauges  and  used  to  calculate  the  axial 
load.  Then  free  vibrations  were  incited  in  the  member  and  were  measured 
by  recording  strain  gauge  outputs  through  the  use  of  an  oscillograph. 

3.1  TEST  MEMBER  AND  SUPPORTS 

The  member  used  in  the  experiments  was  a  mild  steel  1  1/2  x  2  x  1/8 
angle  ninety  inches  long.  The  nominal  thickness  agreed  with  the  actual 
measured  value.  When  the  member  was  clamped  in  the  apparatus  it  was  straight. 

The  different  methods  of  supporting  the  member  for  the  tests  are 
shown  in  figures  3.1,  3.2,  and  3.3.  Figure  3.1  shows  the  experimental 
approximation  of  fixed  supports.  The  member  is  clamped  between  six  inch 
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lengths  of  mild  steel  having  angular  cross  sections  and  equal  legs  one 
half  inch  thick.  For  the  upper  piece  of  the  support  the  angle  has  three 
inch  legs  while  the  lower  piece  has  three  and  one  half  inch  legs.  The 
clamps  were  milled  to  fit  the  member.  The  member  was  placed  in  the  clamps 
as  shown.  Figure  3.2  shows  the  member  bolted  to  the  support  with  the  long 
leg  being  bolted  down.  The  short  leg  of  the  angle  is  shown  bolted  down  in 
figure  3.3.  The  last  two  methods  of  supporting  the  member  are  the  most 
common  for  members  of  angular  cross  section*  but  do  not  approximate  fixed 
ends  as  well  as  the  method  shown  in  figure  3.1. 

3.2  APPARATUS 

The  apparatus  used  in  the  experiments  is  diagrammed  in  detail  in 
figure  3.4. 

An  I-beam  twenty  feet  long,  two  feet  deep  with  nine  inch  wide 
flanges  was  used  as  a  base  which,  compared  with  the  test  member,  could  be 
considered  rigid.  The  apparatus  for  supporting  the  test  member  was  fixed 
to  this  base  as  shown  in  figure  3.4. 

Figure  3.4  shows  the  apparatus  used  for  applying  the  axial  load. 
The  hydraulic  ram  is  mounted  in  the  slide  and  yoke  system.  The  base  of 
the  support  is  welded  onto  a  three  quarter  inch  plate  which  has  its  edges 
milled  to  fit  grooves  milled  out  of  one  and  one  half  inch  bars.  The  bars 
are  welded  to  the  I-beam.  The  grooves  were  greased  to  decrease  friction. 
The  apparatus  was  designed  so  that  the  ram  could  be  mounted  to  put  the 
test  member  into  either  tension  or  compression.  The  volume  of  fluid  in  the 
ram  was  kept  constant  at  any  desired  level  by  shutting  off  a  valve  between 
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the  pump  and  ram.  This  was  done  to  maintain  a  constant  value  of  the 
pressure  and  thus  the  axial  load  while  the  vibration  records  were  being 
taken. 

The  strain  caused  by  the  applied  axial  load  was  measured  by 
eight  Baldwin-Lima-Hami 1  ton  SR-4  strain  gauges.  The  gauges  were  po¬ 
sitioned  three  inches  from  the  support  attached  to  the  slide  when  the 
beam  was  clamped  and  when  the  small  leg  was  bolted  to  the  base  of  the 
clamp.  The  gauges  were  three  inches  from  the  fixed  end  when  the  long 
leg  was  bolted  to  the  base  of  the  clamp.  The  value  of  the  axial  load  was 
calculated  from  the  strain  in  the  gauges.  In  the  calculation  of  the  load 
the  value  of  Young's  modulus  was  assumed  to  be  29.5  x  10  psi.  The  same 
test  specimen  was  used  by  Didrikson  (9)  who  showed  that  the  gauges  used 
were  adequate  to  accurately  measure  the  axial  strain. 

Two  Baldwin-Lima-Hami 1  ton  SR-4  strain  gauges  were  used  to 
measure  the  strain  caused  by  the  vibrations  of  the  member.  The  gauges 
were  placed  forty-five  inches  from  the  end  of  the  test  member  and  at  the 
extremity  of  each  leg.  The  electrical  signals  from  the  gauges  were  ampli¬ 
fied  by  an  Ellis  Associates  BAM-1  Bridge  Amplifier  and  Meter.  The  signal 
from  this  was  recorded  using  a  Brush  Oscillograph,  model  16-2308-00. 

The  axial  strain  was  measured  by  the  BAM-1  which  has  1%  full 
scale  accuracy.  The  frequency  response  of  the  BAM-1  was  0  -  20,000  cycles 
per  second  within  5%.  The  frequency  response  for  the  Brush  Oscillograph 
was  ±5%  for  frequencies  to  1000  cps. 
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3.3  EXPERIMENTAL  PROCEDURE 

For  each  of  the  tests  the  member  was  loaded  in  thousand  pound 
increments  from  zero  to  ten  thousand  pounds  in  tension  and  from  zero  to 
near  the  buckling  load  in  compression. 

First  a  zero  reading  was  taken  for  the  strain  gauges  used  for 
measuring  the  axial  strain.  Then  the  load  was  applied  until  the  gauges 
registered  the  amount  of  strain  caused  by  the  desired  axial  load.  Next 
the  beam  was  struck  in  such  a  way  that  one  of  the  two  lowest  modes  of 
vibration  was  excited  and  recorded.  The  axial  strain  was  then  read  again 
and  the  load  was  released.  Then  the  zero  reading  was  checked  for  drift. 
This  procedure  was  repeated  for  each  load. 


Figure  3.1  BOTH  LEGS  CLAMPED  -  OBLIQUE  VIEW 


Figure  3.2  LONG  LEG  CLAMPED  -  OBLIQUE  VIEW 


Figure  3.3  SHORT  LEG  CLAMPED  -  OBLIQUE  VIEW 
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Figure  3.4  SCHEMATIC  VIEW  OF  EXPERIMENTAL  APPARATUS  (NOT  TO  SCALE) 


CHAPTER  IV 


RESULTS 

This  chapter  presents  the  results  obtained  in  the  investigation 
into  the  effect  of  boundary  conditions  upon  the  coupled  torsional  vibrations 
of  an  axially  loaded  angle  member  with  an  unequal  leg  cross  section.  The 
experimental  results  are  compared  with  the  theoretical  results,  which  were 
obtained  by  programming  the  University's  IBM-360  computer  to  do  the  trial 
and  error  procedure  outlined  in  Chapter  II.  The  required  input  for  the 
program  was  one  value  of  3  and  various  values  of  frequency  squared  and 
thrust.  The  program  results  gave  the  value  of  the  boundary  condition 
determinant  for  a  particular  set  of  input  data.  When  the  sign  of  the  de¬ 
terminant  changed,  the  value  of  the  natural  frequency  was  between  the  two 
values  of  frequency  corresponding  to  the  determinant  values  of  different 
sign.  The  difference  could  be  made  as  small  as  desired  by  reading  the 
values  of  frequency  squared  into  the  program  with  as  small  a  difference  as 
desired.  The  difference  of  the  values  of  frequency  squared  used  for  the 
fixed  ends  was  100  (radians/second)  .  For  the  accuracy  to  which  the  graphs 
could  be  plotted  a  solution  this  close  to  the  actual  value  was  not  required, 
so  in  order  to  save  computing  time  the  remaining  results  were  obtained  with 
a  difference  between  the  two  values  of  frequency  squared  of  1000  (radians/ 
second)^. 

The  experimental  values  of  the  two  lowest  natural  frequencies 
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for  the  values  of  thrust  used  were  easily  obtained  by  the  experimental 
procedures  outlined  in  the  third  chapter.  The  frequencies  higher  than 
the  second  were  difficult  to  obtain  because  of  the  predominance  of  the 
lower  two  frequencies  of  the  member.  No  suitable  method  could  be  found 
to  excite  only  the  higher  frequencies.  For  this  reason  the  theoretical 
results  presented  in  this  chapter  are  the  first  three  frequencies  for  the 
member  with  fixed  ends  and  only  the  first  two  frequencies  for  the  beam 
with  elastic  end  conditions. 

The  frequencies  measured  experimental ly  resulted  from  different 
modes  of  vibration  of  the  member.  This  was  evident  because  they  were 
obtained  by  striking  the  member  from  various  angles  so  that  different 
amounts  of  the  three  displacements  were  present  in  the  vibration. 

4.1  FIXED  ENDS 

The  theoretical  results  obtained  by  using  fixed  end  boundary 
conditions  are  presented  in  table  4.1  and  figures  4.1  and  4.2.  The 
boundary  condition  determinant  used  for  fixed  ends  is  presented  in  Ap¬ 
pendix  1.  These  results  are  presented  to  complete  the  theoretical  re¬ 
sults  of  Didrikson  (9).  The  program  used  in  the  previous  work  did  not 
give  results  for  all  values  of  thrust  for  the  second  and  third  modes. 

The  determinant  solution  was  changed  so  that  the  frequencies  of  vibration 
could  be  obtained  for  any  value  of  thrust  for  the  first  three  modes.  The 
curves  shown  are  those  of  the  first  two  fundamental  modes  of  vibration  on 
figure  4.1  and  the  third  mode  on  figure  4.2.  The  theoretical  buckling 
load  of  the  member  is  the  value  of  thrust  for  which  the  frequency  of  the 
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lowest  natural  mode  is  zero. 

The  results  of  the  test  which  most  closely  approximates  clamped 
ends  are  given  in  table  4.2.  These  results  are  compared  with  the  theoretical 
results  for  the  first  two  fundamental  modes  in  figure  4.3.  This  figure 
shows  that  the  experimental  values  lie  on  a  straight  line  which  is  lower 
than  the  plot  of  the  theoretical  values  and  has  a  smaller  slope.  The 
experimental  line  when  extrapolated  crosses  the  thrust  axis  at  a  value 
lower  than  the  theoretical  buckling  load.  Didrikson  (9)  showed  that  the 
variation  of  the  theoretical  curves  caused  by  varying  the  values  of  the 
modulus  of  elasticity  E  and  the  modulus  of  rigidity  G  would  not  account 
for  the  differences  between  the  theoretical  and  experimental  values  of 
the  frequencies  squared.  When  the  moduli  are  varied  the  plotted  line  of 
frequency  squared  versus  thrust  is  moved  up  or  down  but  the  slope  of  the 
line  does  not  change.  The  differences  between  the  theoretical  and  experi¬ 
mental  curves  are  explained  by  taking  the  end  conditions  to  be  elastic 
instead  of  fixed  ends. 

4.2  ELASTIC  END  CONDITIONS 

The  theoretical  results  discussed  in  this  section  are  those  for 
which  the  boundary  conditions  are  expressed  in  equation  2.31.  For  those 
end  conditions  the  values  of  the  spring  constant  3  were  taken  as  a  para¬ 
meter  and  various  curves  of  frequency  squared  versus  thrust  were  obtained 
for  the  first  two  modes  by  using  the  procedure  outlined  in  chapter  II. 

The  value  of  3  was  varied  from  zero,  to  infinity,  which  corresponds  to 
fixed  supports.  Table  4.3  gives  values  for  and  figure  4.4  shows  representa- 
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tive  curves  of  thrust  versus  frequency  squared  obtained  for  the  various 
values  of  3  for  the  first  mode  of  vibration.  Similarly,  table  4.4  and 
figure  4.5  give  the  theoretical  results  for  the  second  natural  frequency. 
From  figure  4.4  and  4.5  it  can  be  seen  that  the  larger  the  value  of  the 
equivalent  spring  constant  3,  the  greater  the  slope  of  the  curve,  and  the 
higher  the  point  where  the  curve  crosses  the  frequency  squared  axis. 

The  points  shown  in  figures  4.6,  4.7,  and  4.8  show  the  experi¬ 
mental  results  for  the  various  end  conditions.  These  were  used  to  establish: 
effective  values  of  3  for  the  first  two  modes.  The  theoretical  curves 
for  the  values  of  3  that  best  fit  the  experimental  results  are  shown  in 
the  same  figures,  and  can  be  seen  to  fit  the  points  very  well.  The  small 
deviations  are  probably  due  to  errors  in  the  experimental  method. 

Since  the  value  of  the  effective  spring  constant  3  is  smaller 
for  the  first  mode  of  vibration  than  for  the  second  mode  it  appears  that 
the  value  of  3  depends  upon  the  mode  of  vibration. 
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TABLE  4.1 

CALCULATED  NATURAL  FREQUENCIES  OF  ANGLE  UNDER  AXIAL  LOAD 

FIXED  ENDS 


THRUST  FREQUENCY 

FREQUENCY  SQUARED 

1  s  t  2nc* 

3rd 

-jSt  2ncl  3rcl 

lb.  cps  cps 

...  cps 

2  2 

(r/sec)  (r/ sec)  (r/sec) 

-6,000  7.1  39.5  82.5 


-4,000 

17.4 

42.1 

84.5 

-2,000 

23.6 

44.6 

86.8 

0 

28.1 

47.0 

88.4 

2,000 

32.3 

50.8 

90.3 

4,000 

35.8 

51.4 

92.0 

6,000 

39.0 

53.5 

93.6 

8,000 

45.0 

54.8 

94.8 

10,000 

49.4 

57.6 

98.0 

2,000 

61 ,400 

268,000 

12,000 

70,000 

283,800 

22,000 

78,000 

297,700 

31  ,200 

87,200 

308,800 

41  ,200 

95,900 

322,500 

50,700 

104,200 

333,900 

60,000 

113,000 

345,000 

69,200 

118,200 

355,000 

78,200 

131  ,000 

366,000 

000 L  X  (QN033S/SNVICMD  Q3«Vnt)S  A0N3nC)3dd 


Figure  4.1  CALCULATED  FREQUENCY  SQUARED  vs  AXIAL  LOAD  -  FIRST  AND  SECOND  MODE  -  FIXED  ENDS 
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TABLE  4.2 

EXPERIMENTALLY  DETERMINED  NATURAL  FREQUENCIES 
OF  ANGLE  UNDER  AXIAL  LOAD 
CLAMPED  ENDS 


THRUST 

FREQUENCY 
-|St  2nd 

FREQUENCY  SQUARED 

1 s t  2nd 

lb. 

cps 

cps 

2 

(rad/sec) 

(rad/sec) 

-4,000 

15.7 

41 .5 

9,720 

68,100 

-3,000 

19.3 

42.3 

14,700 

70,700 

-2,000 

21.8 

43.5 

18,800 

75,000 

i 

«# 

o 

o 

o 

24.0 

44.7 

22,800 

79,000 

0 

26.6 

45.6 

27,900 

82,000 

1 ,000 

29.0 

46.9 

33,500 

86,500 

2,000 

30.8 

47.9 

37,500 

90,600 

3,000 

31.9 

48.8 

40,200 

93,600 

4,000 

34.0 

50.0 

45,800 

98,600 

5,000 

35.4 

51 .0 

49,500 

103,000 

6,000 

36.8 

51 .8 

53,500 

106,400 

7,000 

38.7 

52.5 

59,300 

109,000 

8,000 

39.9 

53.7 

63,000 

114,200 

9,000 

41.1 

55.6 

66,600 

122,500 

10,000 

42.6 

56.1 

71  ,900 

125,000 

Figure  4.3  PLOT  OF  CALCULATED  FREQUENCY  SQUARED  vs  AXIAL  LOAD  COMPARED  WITH  EXPERIMENTAL  RESULTS 
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TABLE  4.3 

CALCULATED  NATURAL  FREQUENCIES  OF  ANGLE  UNDER 
AXIAL  LOAD  FOR  VARIOUS  VALUES  OF  SUPPORT 
SPRING  CONSTANTS 
FIRST  MODE 


THRUST 

lb. 

6=  0.0 

B=5xl04 

FREQUENCY  SQUARED  (RAD/SEC)2 

3= lxl O5  3=2x1 05  3=4x1 05 

3=1 xl 06 

3  =  00 

-6,000 

-- 

-- 

-- 

-- 

-- 

1 ,000 

2,000 

-4,000 

-- 

-- 

-- 

4,000 

7,000 

10,000 

12,000 

-2,000 

-- 

4,000 

8,000 

12,000 

16,000 

19,000 

22,000 

0 

7,000 

12,000 

17,000 

21 ,000 

25,000 

29,000 

31 ,200 

2,000 

15,000 

20,000 

25,000 

29,000 

34,000 

38,000 

41 ,200 

4,000 

23,000 

29,000 

33,000 

38,000 

43,000 

46,000 

50,700 

6,000 

30,000 

37,000 

41 ,000 

46,000 

51  ,000 

56,000 

60,000 

8,000 

38,000 

45,000 

50,000 

54,000 

60,000 

63,000 

69,200 

10,000 

46,000 

53,000 

58,000 

63,000 

68,000 

73,000 

78,200 
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TABLE  4.4 

CALCULATED  NATURAL  FREQUENCIES  OF  ANGLE  UNDER 
AXIAL  LOAD  FOR  VARIOUS  VALUES  OF  SUPPORT 
SPRING  CONSTANT 
SECOND  MODE 

THRUST  FREQUENCIES  SQUARED  (RAD/SEC)2 


lb. 

6  =  0.0 

3=1x10 

8=2x1 0 

8=4x10 

8=1x10 

8=4x10 

8  =  00 

-6,000 

4,000 

14,000 

23,000 

35,000 

47,000 

57,000 

61 ,000 

-4,000 

12,000 

22,000 

31 ,000 

43,000 

55,000 

65,000 

70,000 

-2,000 

19,000 

31 ,000 

39,000 

51  ,000 

63,000 

74,000 

78,000 

0 

28,000 

39,000 

46,000 

59,000 

72,000 

83,000 

87,000 

2,000 

36,000 

47,000 

57,000 

67,000 

80,000 

91  ,000 

96,000 

4,000 

43,000 

53,000 

63,000 

74,000 

88,000 

100,000 

104,000 

6,000 

51 ,000 

64,000 

70,000 

83,000 

96,000 

108,000 

113,000 

8,000 

59,000 

72,000 

79,000 

91  ,000 

104,000 

119,000 

118,000 

10,000 

67,000 

80,000 

86,000 

99,000 

113,000 

126,000 

131 ,000 

GO 


C\J 


CO 


CO 


•xj- 


CXI 


o 


CXJ 

I 


'^r 

I 


CO 


CO 

o 


X 

-Q 


Q 

< 

o 


X 

c 


a(QN033S/SNViavy)  Q3yVnt)S  A0N3flb3Hd 
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cc 

o 

Q_ 

CL. 

ZD 

GO 


O 

GO 

L±J 

ZD 

< 


GO 

=> 

o 

I— I 

oc 
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cc 
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a 
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X 
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co 

> 

CD 

LU 

cc 

■=c 


ZD 

<zr 

GO 

X 
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TABLE  4.5 

CALCULATED  NATURAL  FREQUENCIES  OF  ANGLE  UNDER  AXIAL  LOAD 
(SUPPORT  SPRING  CONSTANTS  CHOSEN  TO  GIVE  BEST  FIT  OF 
CALCULATED  FREQUENCIES  TO  EXPERIMENTAL  FREQUENCIES) 

CLAMPED  ENDS 


THRUST 

FREQUENCY  SQUARED 
FIRST  MODE 

(RAD/SEC)2 
SECOND  MODE 

lb. 

3  =  8  x  105 

3  =  4  x  106 

-6,000 

500 

__ 

-4,000 

9,000 

65,000 

-2,000 

19,000 

74,000 

0 

27,000 

82,000 

2,000 

37,000 

91  ,000 

4,000 

46,000 

100,000 

6,000 

55,000 

108,000 

8,000 

64,000 

117,000 

10,000 

73,00 

125,000 
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TABLE  4.6 

EXPERIMENTALLY  DETERMINED  FREQUENCIES  OF  ANGLE 

UNDER  AXIAL  LOAD 
LONG  LEG  CLAMPED 


THRUST 

FREQUENCY 

1st  MODE  2nd  MODE 

FREQUENCY  SQUARED 

1st  MODE  2nd  MODE 

lb. 

cps 

cps 

2 

(rad/sec) 

2 

(rad/sec) 

-4,000 

15.8 

37.5 

9,860 

55,500 

-3,000 

18,7 

39.6 

13,800 

61 ,900 

-2,000 

20.5 

41.1 

16,600 

66,700 

o 

o 

o 

0% 

1 

23.2 

42.3 

21 ,300 

70,600 

0 

24.8 

42.9 

24,250 

72,600 

1  ,000 

26.9 

44.1 

28,500 

76,600 

2,000 

30.3 

45.4 

36,200 

81 ,600 

3,000 

31.4 

46.4 

38,900 

84,900 

4,000 

33.2 

47.6 

43,500 

89,400 

5,000 

35.0 

49.3 

48,300 

95,800 

6,000 

36.1 

50.0 

51 ,400 

98,600 

7,000 

37.7 

51.4 

56,000 

104,000 

8,000 

38.6 

52.0 

58,700 

106,500 

9,000 

40.3 

53.2 

64,000 

111  ,600 

10,000 

41  .9 

56.8 

69,300 

127,000 

56 


TABLE  4.7 

CALCULATED  NATURAL  FREQUENCIES  OF  ANGLE  UNDER  AXIAL  LOAD 
(SUPPORT  SPRING  CONSTANTS  CHOSEN  TO  GIVE  BEST  FIT  OF 
CALCULATED  FREQUENCIES  TO  EXPERIMENTAL  FREQUENCIES) 

LONG  LEG  CLAMPED 


THRUST 

FREQUENCY  SQUARED 
FIRST  MODE 

(RAD/SEC)2 
SECOND  MODE 

lb. 

6  =  4  x  105 

3  =  1.2  x  106 

-6,000 

-- 

% 

49,000 

-4,000 

8,000 

57,000 

-2,000 

16,000 

66,000 

0 

25,000 

74,000 

2,000 

34,000 

82,000 

4,000 

43,000 

91 ,000 

6,000 

52,000 

99,000 

8,000 

60,000 

107,000 

10,000 

68,000 

116,000 

* 
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TABLE  4.8 

EXPERIMENTALLY  DETERMINED  NATURAL  FREQUENCIES  OF 
ANGLE  UNDER  AXIAL  LOAD 
CLAMPED  ENDS 


THRUST 

FREQUENCY 

1st  MODE  2nd  MODE 

FREQUENCY  SQUARED 

1st  MODE  2nd  MODE 

lb. 

cps 

cps 

2 

(rad/sec) 

2 

(rad/sec) 

-4,000 

17.0 

37.2 

11,450 

53,700 

-3,000 

19.6 

39.1 

15,100 

60,600 

-2,000 

21 .4 

40.5 

18,000 

64,800 

-1 ,000 

23.5 

41  .8 

21 ,550 

69,000 

0 

25.7 

42.8 

26,100 

72,400 

1 ,000 

27.7 

44.1 

30,300 

76,700 

2,000 

31 .0 

45.3 

38,000 

80,900 

3,000 

32.0 

46.4 

40,600 

84,600 

4,000 

34.0 

47.3 

45,800 

88,200 

5,000 

35.0 

48.3 

48,400 

91  ,900 

6,000 

36.5 

49.6 

53,800 

97,400 

7,000 

38.1 

50.7 

57,100 

101,000 

8,000 

39.3 

51  .9 

61  ,000 

106,200 

9,000 

40.4 

52.7 

64,500 

109,600 

10,000 

41.7 

53.6 

68,600 

113,500 

59 


TABLE  4.9 

CALCULATED  FREQUENCIES  OF  ANGLE  UNDER  AXIAL  LOAD 
SUPPORT  SPRING  CONSTANTS  CHOSEN  TO  GIVE  BEST  FIT  OF 
CALCULATED  FREQUENCIES  TO  EXPERIMENTAL  FREQUENCIES 

SHORT  LEG  CLAMPED 


THRUST 

FREQUENCY  SQUARED 
FIRST  MODE 

p 

(r/sec) 

SECOND  MODE 

lb. 

3  =  6  x  105 

3  =  1.1  x  1 06 

-6,000 

-- 

-- 

-4,000 

8,000 

56,000 

-2,000 

18,000 

64,000 

0 

27,000 

73,000 

2,000 

36,000 

81  ,000 

4,000 

44,000 

89,000 

6,000 

53,000 

98,000 

8,000 

62,000 

106,000 

10,000 

71 ,000 

114,000 

I 

CHAPTER  V 


CONCLUSIONS  AND  APPLICATIONS 


5.1  CONCLUSIONS 

The  trial  and  error  procedure  of  Gere  and  Lin  (8)  used  to  solve 
the  three  simultaneous  partial  differential  equations  describing  the  vi¬ 
brations  of  a  thin  walled  member  of  open  cross  section  will  work  for  any 
set  of  elastic  boundary  conditions. 

When  the  boundary  conditions  of  a  member  are  known  and  are  ex¬ 
pressed  in  terms  of  the  elastic  spring  constants  of  Chapter  II  then  the 
axial  load  on  the  single  span  member  can  be  predicted  by  measuring  the 
lowest  or  the  second  lowest  natural  frequency  of  vibration. 

5.2  POSSIBLE  APPLICATION 

This  investigation  of  the  effect  of  end  conditions  upon  the 
frequency  of  an  axially  loaded  beam  was  carried  out  in  an  effort  to  find 
a  method  of  determining  the  axial  load  in  a  member  of  an  existing  structure 
The  results  presented  show  that  if  the  end  conditions  of  a  member  are  known 
then  the  thrust  on  the  member  can  be  determined  accurately  by  measuring 
either  the  first  or  second  frequency  and  reading  the  value  of  axial  load 
from  the  theoretical  curve  of  frequency  squared  versus  thrust.  This  would 
give  an  experimental  method  for  confirming  the  value  of  load  used  in  design 
ing  a  member  in  a  complicated  structure. 

This  application  depends  upon  knowledge  of  the  end  conditions. 
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and  in  the  case  above  one  would  need  to  know  the  value  of  3  for  the  lowest 
frequency.  The  value  of  the  3  was  seen  to  depend  upon  the  mode  of  vibra¬ 
tion,  so  the  values  of  the  spring  constant  would  have  to  be  found  experi¬ 
mentally  for  different  end  conditions.  An  attempt  was  made  to  find  a 
value  of  3  statically,  but  the  result  was  not  the  same  as  either  of  the 
values  of  3  for  the  first  two  modes.  Thus  in  order  for  the  method  out¬ 
lined  in  this  thesis  to  be  of  practical  importance,  it  will  be  necessary 
to  determine  the  relationship  of  the  spring  constants  to  the  type  of  end 
support. 
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APPENDIX  I 


In  chapter  II  a  solution  was  presented  for  the  fundamental 
modes  of  vibration  of  an  axially  loaded  member  of  general  cross  section. 
This  was  reduced  to  the  specific  case  of  an  angle  member  with  unequal 
leg  cross  section.  The  general  case  considered  in  Chapter  II  reduces 
to  the  case  of  a  member  with  simple  supports  or  clamped  supports  de¬ 
pending  upon  the  value  of  the  spring  constant  3.  In  this  appendix  the 
boundary  condition  determinant  will  be  derived  for  the  case  of  fixed 
ends  and  the  problem  will  be  solved  without  the  use  of  the  computer  for 
simple  supports.  Both  of  the  above  solutions  were  used  to  check  the  re¬ 
sults  from  the  computer. 

CLAMPED  SUPPORTS 

The  method  used  to  solve  the  problem  with  clamped  supports  is 
the  same  as  the  solution  for  general  supports  in  chapter  II,  save  that 
the  proper  end  conditions  had  to  be  employed.  This  makes  the  boundary 
condition  determinant  slightly  different.  The  boundary  conditions  for 
the  member  with  both  ends  clamped  are 


v  =  w  =  <p  =  0 

dv  _  dw  _  d(j)  _  n 
dx  ”  dx  dx 

This  follows  from  the  discussi 


For  x  =  0  x  =  i 
For  x  =  0  x  =  Z 

n  of  boundary  conditions  in  the 
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■ 
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second  chapter.  It  has  also  been  shown  that  the  solutions  of  V,  W,  and 
$  may  be  expressed  in  the  form 


rlx  "rlx  1  1 

V  =  Ai e  +  A2e  +  cosr^x  +  A^  sinr^x  + 


rlx  -r-|X  ,  1 

W  =  y-| A-j e  +  y-jA^e  +  y2A3  cos +  y2A4  sinr^x  +  . 

rlx  ~rlx  1  1 

$  =  \p.| A-| e  +  ^-j^2e  +  ^2^3  cosr2x  +  ^2^4  s''nr2x  +  •••• 

Now  to  obtain  the  boundary  condition  determinant  for  fixed  ends, 
these  values  of  V,  W,  and  $  are  substituted  into  the  boundary  conditions 
above  and  the  determinant  of  the  coefficient  is  set  equal  to  zero.  This 
determinant  was  earlier  called  the  boundary  condition  determinant.  For 


clamped  ends 

it  is  of  the 

form 

1 

1 

1 

0 

r,£ 

e 

-r-|£ 

e 

cosr2£ 

sinr2& 

Y1 

Y2 

0 

/l* 

Y-|e 

-r-jH 

Y-|e 

Y2cosr2£ 

Y2sinr2£ 

*1 

*1 

\p2 

0 

r,H 

-r,& 

e 

ip^cos  r^i 

^2sinr2£ 

ri 

"rl 

0 

r2 

r-|£ 

r-je 

-r-j  9j 

-r-je 

0 

r2 

Vi 

"Ylrl 

0 

Y2r2 

r,£ 

W 

-r,  & 

-vie 

-Y2r2sinr2& 

Y2>^2cosr2£ 

Vi 

'^iri 

0 

^2r2 

r-|£ 

^irie 

-r^ 

-^irie 

-^2r2s 1 nr2^ 

^2r2cosr2£ 

=  0 
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To  solve  for  the  natural  frequency  of  a  member  with  clamped  ends 
this  determinant  is  used  in  the  trial  and  error  procedure  outlined  in 
chapter  II.  In  the  case  where  the  warping  constant  cw  is  zero  the  eleventh 
and  twelfth  rows  and  columns  of  the  above  determinant  are  deleted. 

SIMPLE  SUPPORTS 

The  natural  frequencies  of  a  member  with  simple  supports  could 
be  solved  for  by  the  method  outlined,  but  they  can  also  be  obtained  by 
hand  very  easily  as  was  done  by  Vlasov  (2).  The  following  shows  the 
method  used  to  solve  the  simple  support  case  by  hand.  The  boundary 
conditions  are 


v  =  w  =  $  =  0 


2  2  2 

d  v  _  d^W  _  dS 

2  2  2 

dx^  dx  dx^ 


=  0 


These  support  conditions  suggest  the  use  of  the  solutions 


w  w  •  nirx 
V  =  V0  sin  — 


..  .  nux 

W  =  W0  sin  — 


,  •  nnx 

4.  =  sin  — 


Substitute  these  values  for  V,  W,  and  $  into  equations  2.22, 
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2.23,  and  2.24  and  simplify  the  results  by  observing  that  the 
natural  frequencies  for  a  simply  supported  member  are  given  by 


,  2  _  /mM 
y  mA  v  ; 


2  _  EIn  /nTTs4 
z  mA  '  i  ' 


P 


2 

<P 


nTT\2 


<x> 


2  2  2 
n  7T  Ec  +  £  Gc 

( - ^ - ) 


£2ml 


0 


Further  simplification  is  made  by  putting 


k 


2 

n 


P  /Httn  2 
mA  ^  i  ' 


From  this  it  is  seen  that 


P 


n 


P  /n7T\2 
mA  '  i  ‘ 


The  above  substitutions  and  simplifications  reduce  the  solution 
of  the  differential  equations  to  the  solution  of  the  following  determinant 


2  b  2 

Py  -  kn 
0 


A  ck2 
!0  z  n 


0 


2  is  2 

P  -  k 
vz  n 


fck  2 
J0  y  n 


c  k 
z  n 


-ck' 
y  n 

2  ,  2 
p ,  -  k 

n 


=  0 


•  . 


. 
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which  gives  the  following  cubic  equation 


unV-^[py2  n-^]  +  pz2  [i  -  +  p>n2)2 


+  [P  2P  2  +  P  2P  2  +  P  2  P  2]  k  2  +  ~  P  2P  2P  2  =  0 
Ip  z  y  4>  y  n  I  $  y  z 


? 

This  equation  was  solved  by  hand  for  kn  and  thus  for  the  squares 

2 

of  the  natural  frequencies  pn  as  a  check  upon  the  general  solution  ob¬ 
tained  in  chapter  II. 


. 


APPENDIX  II 


The  method  of  solution  outlined  in  chapter  II  was  programmed 

for  a  computer  in  the  Fortran  IV  language.  Two  programs  were  used  to 

obtain  the  natural  frequencies  of  the  member  for  various  support  conditions 

with  an  applied  axial  load.  The  first  program  uses  the  member  constants, 

2 

value  of  thrust  P,  and  frequency  squared  pn  to  set  up  and  solve  the  fifth 
degree  characteristic  equation.  The  output  is  the  square  of  the  roots  for 
the  equation.  The  second  program  uses  the  member  constants,  values  of 
thrust,  frequency  squared,  spring  constant  3,  and  roots  of  the  characteristic 
equation  to  establish  and  evaluate  the  boundary  condition  determinant.  The 
program  writes  the  value  of  3,  thrust,  frequency  squared,  and  of  the  de¬ 
terminant  as  output. 

To  obtain  the  natural  frequency  the  value  of  the  determinant  is 
plotted  versus  frequency  squared  for  particular  values  of  thrust  and  spring 
constant  3.  The  natural  frequency  squared  is  where  the  line  crosses  the 
frequency  squared  axis. 

The  two  programs  used  are  listed  in  this  appendix.  When  applying 
these  programs  to  other  problems  it  should  be  noted  that  they  are  for  an 
unequal  leg  angle  cross  section,  which  means  the  value  of  cw  is  zero.  For 
a  general  cross  section  cw  is  not  zero  and  the  programs  must  be  changed. 

The  solution  for  the  roots  is  only  for  a  fifth  order  equation. 

The  following  is  a  list  of  notations  used  in  the  programs 
HA  A,  area  of  the  cross  section  square  inches 

HC  c,  torsion  constant  in  (inches)4 
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. 


polar  moment  of  inertia  about  shear  center  in  (inches) 

2 

1  h  qpr 

mass  density  in  — a — 

length  of  member  in  inches 
eccentricity ,  in  inches 

principal  centroidal  moment  of  inertia  about  n  axis 

/  x  4 

in  (inches) 

modulus  of  elasticity 

2 

square  of  circular  frequency  in  (rad. /sec) 
modulus  of  rigidity  in  psi 
thrust  in  pounds 

polar  moment  of  inertia  about  centroidal  moment  of 

/  v  4 

inertia  about  the  l,  axis  in  (inches) 
y  coordinate  of  centroid  in  inches 
z  coordinate  of  centroid  in  inches 


; 

J 

PROGRAM  ] 


o  T  MP  N  F  T  ON  A  (  *  0  )  ,  P  (  ^  o  )  ,r(tn)  ,  R  P  (  R  o  ) 
pr  An  (  F  ,  2  )  HA  ,  f  I  9  ZI  » nr  T  *  9>I  •  Ff  7  .  ROY  ,Hr  ,fy  , 
2  FORMAT  (  1  X  *7F10.5/1  X  *F9.'5  *3X  *2F12 .2  »F9.6 
wp  TTF(G,4)HA,FI  *Z  I  ->HOT  ,  F  I  *FCZ  *FCY  *  HC  »  r  Y 
4  FORMAT  (  1  X  ,7 FI  0.5  *  /]  X  ,4F1F.F  *F1;n.  2  ) 

V A  =  F  Y # R Y * Z  I  * F  I  /  (  DIS|*D-N*DN*HA*H4*HC  I  ) 
VB=G*HC 

Vf  =F Y / ( DN*ON*HA*HA ) 

\/n  =  c  I  /HA  -  (  F  T  *FC7*PG7+7  T  *FCY-*FCY  )  /HO  T 
VE=EY*Z I *E I*S I /HOI 
VG  = 1 . / ( DN*DN*HA ) 

100  RFAn(F,7P)F,F 
'VR  ITF(F%?F)P»F 
PR  FORMAT ( 1X*2F10. 1 > 

C  F  A  =  V  A  * ( S*S I /HA-VR ) 

v  I  =f*o/QN*VD“VP*F/BN 

VJ  =  S*.F*  (  F/HA-VR/HCI  )  /  (  DN*HA  ) 

VK  =p Y* ( 2 • *S*VD— VR ) 

VL=5/DN/HA* ( 2  .*S/HA-2 .*VR/HC I > 

VM=  (  -a  .*F/HA~VR/HOT  ) /ON  # 

A ( 1  ) =1  .  n 

A  (  2  ) =V0* ( V  I -P* VF ) /CF  A 
a  (  o  )  =v.G*  (  VJ-F*VK  )  /  OF  A 
A  (  4  )  =  —  F  /  D  N  *  (_VL-F*FY*VD  )  /CF  A 
A  (  2 ) =  F * F * V M /CF A 
A ( A ) =-F*F*F/CFA 

Nj  -  K 

FP^FL.  =  +  0 . 1  OOP-1  o 

N  =  M+1 

P  =  o .  5 
<=;  p  =  o  ,  o 

G  =  o ,  o 

6  D07I=?  *NJ  tl 
R()  )  = A  (  1  ) 

d  (  7  )  =  A  (  p  )  -  P  *  P  (  1  ) 

7  R ( I ) sA ( I ) -P*p ( T “I ) -G*P ( T -P ) 

|_  =  N  -  2 

n op  T  =  2  * L  9 1 
0(1  )  =  R (  1  ) 

C ( 7 ) =R ( ? ) -P*C ( 1 > 
p  r (  I  )  =  R (  I  ) -D*r (  I - 1  ) -G*c (  1 -0 ) 
c  (  n1  —  i  )  =  o«o  —  P*C  (  N  —  P  )  —  G  *  C  (  N  —  0  ) 
n=r { N - 2 ) **2-C ( N-l ) *C ( N~p 1 
I F ( n . FO . o . o ) GO  TO  1° 
so  OF! P= (P ( N-1 )*C (N-2 )-n (N )*c (N-2 ) ) /o 
OF LG=  (  R  (  N-l  )  *C  {  N-l  )  -p  (  N  )  *C  (  N-2  )  )-/r> 
P=P+DFLP 
G=G  —  DFLG . 

|p  (  A, D  F  (  P F L G  )  •  G T> •  p P  F F  L  )  G 0  TO  A, 
o  a  n  -  p  *  *  o  -  4  .  o  *  G 

TF ( AP.LT .6.0 ) GO  TO  IF 


G  ,RN  ,  OF 
*F9.1  ) 

,G,OM, OF 


. 

' 


in  p p  =  a p * * n , 5 

YY = ( p „ n-P+RR ) /? . n 
YY= ( n . n-P-RP ) /?  .  n 
R  P  (  N  )  =  X  X 
RR ( N-l ) = Y Y 

'•'R  ITF  (  6  *  1  1  )  XX  9  YY  • 

11  FORMAT ( IX* 7H  ROOTS  ?  ?  F 13  .  S ) 
i  l  n  m  s  |s)  —  y 

IF(N*LT«3) GO  TO  17 

i  ?  b  n  )  =  l .  o 

R ( ? ) =A ( ? )-P 
A ( ? ) =R ( ?  ) 

no  i  =  n , n  ,  i 

r  (  I  )  =  A  (  I  )-P*R(  T  - 1  )-G*:R(  I  -7  } 
17  A  (  I  )  =R  (  I  ) 

IF  (  M  .GT  .  3  )  GO  TO  3 
14  P= A ( ? ) 

G=  A ( 3  ) 

GO  TO  Q 
1  5  IT  —  A RS  (  AR  ) 

RR=SQRT ( T ) 

IF  (':B  R  •  G  F  •  •  1  0  F  - 1 0  )'$0  TO  3? 

31  ,'-'R  ITF  (  6  »37  ) 

7q  IFORMAT  (  1  X  » 0  H  TMAG  RTF 
GO  TO  100 

a?  roMP=  (  n  .  o-P  ) /?,.  o 

1  3n  RR  (|S)  )  rfOMP 

■••'p  I  TF  (  6  » 1  6  )  COVP 

16  FORMAT  (  1  X  ? 7 H  D.ROOT  ,E'18.8> 
WR  I T  F ( 7  ? 1 i 1  )  ( R  R (  I  )  *  I  -  7  » 6 ) 

.111  FORMAT  (  IX  »3E15.8  ) 

GO  TO  110 

17  IF  (  N  —  ?  )  1.00>1P?100 

i p  roMP=o . o-A ( ? ) +p 

GO  TO  ]30 
IP  p-p+p 
'G  =  0  +  R 

TI  F  (  P  *  LT  .  1.00  „  )  GO  TO  6 
?n  MR ITF ( 6  « 7  1  ) 

?’l  FORMAT  (  3  X  *  1  3  HP  IF  OV/FR  100 
R  =  —  0 . 3 
P  =  —  O  „  3 
G  =  —  o  0  3 
GO  TO  6 
FND 


PROGRAM  ? 


OOl  1  R  |  F 

opr r  j  c  j  o A) 

A  A  (  2  0 

OOIIPLF 

RRFG  T  FT  ON! 

PT  (K  ) 

OOt  IPl  F 

p  r  f  c  T  S  T  O  M 

0(11, 

n  o  u  p  l  E 

p  R  F  C  T  g  I  0  N 

H  A  ,  F  J 

OOl  | R L F 

RRFG  T  R ION 

A  G  ,  A  P 

OOl  JR|_F 

p  p  F  G  T  c  T  0  M 

7(10) 

R  f  a  n  (  F 

,  2  1  H  A  ,  F  J  , 

7  t ,hg t 

,  n 


M I  1 1  TOP 


11  )  ,XM(  1  1  )  ,nrT(  1  1  )  ,r>'FT.A  (1.0  ) 

,7  T  ,Hf  T  .AT  , rrr  7  ,  FGY  ,HG,frY,  G  ,  RN 

,  AO  ,  A  JK  .  TFMP  .  VI.  A  R'G  ,dFG  T  R 


2  FORMAT ( 1 X  » 7  F 1 0,6/1 X 


CO 


T  ,  rrr  7  ,  F G y  .  t- 1  f  ,  r  Y  ,  G  ,  F>M  »  nF 
,7X  *7F1  2.?  *F9«.6  ,P0.1 


M ,  =r  1  0 

N  M  —  4 

AR=FY*Z I /DM/HA 
A  D  =  F  Y*F I / RN /H A 

IPO  R  F  A  D  (  A  ,  4  )  s  ,  F  9-  R  T  (  1  )  .  R  T  (  2  )  ,  R  T  ( 2  )  ,  R  T  (  4  )  »  7  T  (  6  ) 

A  FORMAT  (  ")  X  9?F  1  0*1/1  X  <CF1  F  „  2  ) 

T  J  - 1 

R  R  T  a  f  TJ)  rGAnpono  +  R  A 
Z ( T J ) =1 . /PFTA ( I J ) 

UIR  T  TF  {  6  4  6A7  1  F  FT  a  (  T  J  ) 

667  format ( 1 X  * D1 5 . R ) 

A  r  =  r  /  R  M  /  H  A 
no  2F  1  =  1%  F 

TOR  f  J  )  =  A  R-X-RT  (  T  )  *RT  (  T  )  +  A 0 * R T  (  T  )  “r 

p>  F  N  (  T  )  =  F  r  7  *  (  A  C  *  R  t  (  I  )  -  F  ) 

TF  {  OFM  (  T  j  *FQ  ,  n  ,  n  ]  r’0  Tn  200 
in  a  o  h  I  (  I  ) =  T  OP (  T  )  /;%'F  M  (  I  ) 

o A M  (  T  ) =  A  0* R  T (  J  )  *P T  (  T  )  +AO*RT  (  T  )  -c 

JR  (  o  am  (  I  )  .  FO  ,  0  •  0  )  GO  TO  ^02 

TT  =  TAM(  2  2 • 6  *  7  *  1 41  F R2 /I  PR,  ) 

12  G AM (  T  ) =  TOP (  J  ) *FGY /  ( G AM (  T  ) *Ff 7  ) 

T  F  (  R  T  (  T  )  )  1  F  ,  1  4  %  1  A 

U  OF  T  (  N  )  =  6  6  6  •  6 

GO  TO  of 

IF  PT (  T  )=-RT (  !  ) 

R  r  (  T  )  =R  T  f  T  )  **  *  F 
R (  I  )  =RT  (  T  ) 

P  (  1  ,  2  *  T  —  1  )  =  1  . 

R  (  1  *2*  I  )  =0* 

r  (  2  ,  o  *  T  —  1  ) =OGOS ( Q  T ( T  )  *  n  G  ) 

D  (  7  ,  2  *  1  )  =  R)  R  T  M  (  R.  T  (  T  1*001 
Ff  l  2  ,  2  *  T  -  1  )  =  G  A  M  (  I  ) 
p  (  2  ,2*1  )  =  O  . 

P  (  4 , 2*1 -1  ) =  G A M (  I  )  * R ( 2 , 2 *  T - 1  ) 
n  •(  4 , 2*  J  )  =G  A  M  (  T  )  *R  (  2  ,  ?#  T  ) 

R  (  F  , 2*1-1 ) =ApHI ( I ) 

P  ( F  ,2*1  )  =  O  * 

d  (  A  .  2*  T - 1  )  =  A  PH !  (T  1  *  R (  2 , 2  *  I “ 1  ) 
d  (  ft  ,  2  *  T  )  =  A  R  H  T  (  T  }  *n  (  2  ,  2  *  T  ) 
d(7, 2  *  T  —  1  )  = - r  Y  *  7 T*PT  (  T  ) **2 
p  (  7 , 2*1  )t-RT(  T  1  * n  F  T  A  (  T J  ) 


R(6) 
n  c 


* 

d(  q,o*t-1  )  —  F  Y  *  7  T  *  P  T  (  T  )  -R  T  (  I  )  *  n  )**FTA  (  T  J) 

n  (  0  ,p.X-T  )  a_FY*ZI*«T  (  T  )  **?*p  (7  ,-?*T  )+PT  (  T  )*«  (  9  ,9*1-1  i^PPTA  (  I  J  1 

r ( q , - 1  ) =~F y*F I *G A M (  I ) *RT ( I  )  ** ? 

p  (  o  ,  9*  T  )  =-PT  (  T  )  *r,AM(  T  )  -SPRT  M  T  ,J  ) 

P(10,«1-1)  -  —  RAM  (  T  >*(  pY*FI*(  RT(  T  )**?  1*P  (  ?  >+RT  (  T  >  WTM  T  J 

?p  (  T  J  )  =-n& v  (  T  )  *  (  FV*FI*  |.pT  (  T  )  Ma'iW1!  9  ,7*  T  )-PT  (  T  )  *Dr"T  A  (  T  J  '  '  * 

9  P  (  9  ,  9  *  T  -  1  )  ) 

(RO  TO  9R 

16  RT ( I ) =PT ( T ) **. * 

’R  (  I  )  =RT  (  I  ) 

P  (  1  ,7*1-1  ) = ]  . 
p  (  1  ,  7  *  T  )  =  1  . 
rvR=np*RT  (  I  ) 

T  p  (  P  r  .  0  T  *  1  7  0  )  GO  T  0  a,  n  7 
p  (  9 , 7*  j  _  i  )  =n':Xp  (  pp*RT  (  I  )  ) 
r(7,7*I)=1./R(?>?*T-1 ) 
n ( 9  ,  7*t -1  ) =  GAM  (  I  ) 
p  (  o  ,  7  *  I  )  =  G  A  M  (  I  ) 

R' (A, 7*1  —  1  )  ~Ci  A  M,  (  I  )  *  P  (  ?  ,  7  *  T  —  1  -1 
P  (  4 , 7  *  I  )=GAM(Ij)*R(?»  ?  *  T  ) ; 
p  (  F  .  9  *  T  —  1  )  =  A  P  H.T  (  T  ) 

P  (  F  ,  7  *  1  )  =  A  P  H  I  (  I  ) 
n  (  A  )  =  A  D'H  T  (  T  1  *P  (  9  ,R*T-1  ) 

p  (  A  ,  7.*:T  )  =  A  P  M  l  (  T  )  *R  (  9  ,  9*  I  ) 

p  (  7 , 9*1-1  )  =PY*7  I  *  (  PT  (  T  )  **?  )  -DT  (  T  )  * D  F T  A*'f  I  J  ) 

p  (  7 ,7*T  )  -FY*7.  T  *  (  RT  (  I  >  **7  >+RT  (  T  )  *RFGT  A  (  KJh 

?*T.1  )  =  n  (  9 , 9*T  -  i  )*  (py*7  Tf-RT  (  T  V‘**7+RT  (  T' )  *p  R  T  A  (  I  J  )  J 
9*T  )  =P  (  9 , 9*  t  )  *  (  FY*7  T  *RT  (  T  )  *'-*7  — RT-(  T  )  ^'pp  l  A-(  T  J  )  > 

P*T_1  )  =PAM  (  T  )  *  F  Y  *  P  T  ,*  (  R  T  (  T  1**7  )  -RT  (»IO*GA^(  T  )  *  n  F  T  A  (  T  J  ) 

9*T  )=OAM(  T  )  *  (  p  v  *  r  J  *  P  T  (  T  )  **9+-PT-(  T  1  *  i">  F  T  A  (  T.J  )  1 

)-(rY*q^RT(  I  )**9+RT(  T  )  *  Q  F  T  A  (  TJ)  )*RAM(  T  J*p(7«9*T 


9  F 
Ann 

9  00 

9  0  1 


P  (  o 
n  (  o 
d  (  o 
n  (  O 

pI^’,9*!  )  =  (FY*;FT*RT(  T  )  *  *  7  R  T  (  T  )*^fTA  (  T  J  V )  *G  AM  {  T  )*o 

rOMT  TNUF 
00  7  no  T  ss  1  *  N 
HO  700  J  =  1  M 
a  A  (  T  ,  J  )  = p  (  T  « J  ) 

FALL  0  F  T  F  R  m  (  A  A  «  N  •  0  ) 

WP  J  T  F  (  6  »  ?  0 1  )  D 

format  (AX.OAH  SMPPOIITTOR  AMR.  >  n  1  *  •  B  > 

N  7  =  M  -  7  . 

ooqO  K=1  ?N2 


■1  ) 


9 , 9*T ) 


Y  1  -K+l 
Y.  ?=K'  +  ? 

I  A  P  G  -  ^  1 

\/|  ARG  =  R  (  K  ,K  1  ) 

\  / 1  a  P  O  =  A  P  R  (  V  L  A  R  G  ) 
no A 9  J=K?*N 
AJk'=P(J»^  ^ 

T  F ( VL ARG .GF - 0 APS ( A  JK )  ) GO  TO  6° 
A A  VL ARG  =  DARS ( A JK  ) 

L  A  P  G  =  J 


• 

' 

68  CONTTMtlF 

IF  (  L  A  R  G  —  K  1  )7o,?o.f,7 
6  7  00  6i8  I  =  K  ?  N 

T  p  M  P  =  R  (  T  » K I  ) 

R (  I  , K 1  ) —  P ( T *LARG ) 

6 p  R (  T  ,  LAPG ) bTFMP 
7 r  TF(R(X,i<'1).r"0.,R.O)  GO  TG  7  0  4 
71  RFC  TP  =  1  .  / R ( K  ?  K 1  ) 

0077  j  =  i<  7  ,\l 

77  XM { J ) =-R CK 9 J )*RFG T P 
nopo  I=K 
HOPo  j  =  X  P  *fSJ 

p,  0  r  {  j  ,  J  )  =  R  (  I  *J)+XM(J)*R(  I  » K 1  I 
T  F ( L.ARG-K1  )  p  f  .  P  F  9  Q  1 
81  HOP?  J=1  ♦  N 
T  F  M  P  =  P  (  K  1  »  J  ) 
p ( X  1  9 J ) =R ( L  ARG  9 J  ) 

P p  R  (  I..ARG  ,  J  ) -TFMP 
pc  noon  j  =  1 » N 
p  ('  |  M  =  0  • 

0088  L-K?9N 

p  p  ci  )M  =  FUM  +  XM  ( l.  )  *R  (  L  9  J  ) 

^0  r ( y 1 , J ) *p ( K1  9 J ) -SUM 
n  F  T  (  1  )  =  R  (  1  «1  ) 

00170  m  =  ?  9 M 
0  F  T  m  =  n  (  m  ,  1  ) 

00  110  I  *=  7  9  M 

Tin  hftm  =  R ( v  ,  I  ) *OFT (  T- 1  I -OFTH*R (  I “1  *  1  ) 

17  0  r>FT(VI)#DFTM 

oc  ui  R  J  t  F  (  6  9  Q  6  )  c  t  F  ,  0  F  T  (  M  ) 

on  FORMAT! 1 X  9  7 F 1 0 . 1 9  F 1 F  »  F I 

no  jn  ] no 

PRO  WR  ITF ( A , ) F  9  F 

p  n  1  f  O  R  M  /\  T  (  1  X  9  0  H  0  F  N  (  T  )  =  n  «  7 1  0  •  0  1 
GO  TO  100 

P  0  p  W  R  T  T  F  (  6  9  7  o  p  )  c  «  F 

pro  fORmAT(1X9°H  G a M (  T  ) =  O  ,?F1 O,  1 
GO  TO  100 

P04  v.ir  TTF  (  6  9  Q0F  )  p  9  F 

Prc  rrpmaT(1X»1GH  R  K  *K1  ) =0  9  7  c i  o . i  ) 

GO  TO  100 

nop  \,i  p  i  t  F  (  6  9  6  °  8  )  0  S  »  P  T  (  I  )  9  R  p 
n rp  FORMAT! 1X9001 f*8) 

GO  TO  100 
FM  O 

PUP ROUT  I  ME  OF  T  F  R  M ( A  A  9  N  9  D ' 

nOUPLF  PRFG TFTON  A  A (  70 9?0 )  9  A  ( 70  ,  7n ) 

00  MS  L  F  pDFCT  FION  R  9  ■  9 r  9 ' 

no  7  0 0  I  =  1  9  N 

no  7 gg  J  =  1 9 N 

2  0  0  A  (  I  «  J )  =  A  A (  I  9  J ) 

0=1  9 


1  <K=K  +  1 
T  S  =  K 
TT=K 

p  =  HAPS ( A (K  >K  )  ) 
no  7  f=K,N 
no  7  j  =  y  ,  m 

IF  (  n  A  P  P  (  A.  (  T  ,  J  )  )  ,  l_  F  .  R  )  r-  0  TO  7 
I  c-  =  T 
T  T  =  J 

o  =  n  A  PS ( A  (  I  »  J  )  ) 

7  OOMT  T  Mi l(r 

TF  (  I  c.LF.m  00  T o  ? 
no  77  J  =  K  *  N 
r = a  ( i  s  *  J  ) 

A  (  I  S  »  J  )  =  A  (  K  *  J  ) 

^7  A ( <  ,  J ) =-C 
7  TF  (  IT.lrF.K )  00  TO  4 
no  4  7  T  =  K  •  M 
r= a ( i , r  t  ) 

A  {  t*IT)=A(T,K) 

4 7  A (  I  »  K ) =-r 
4  n=A(K*K)*D 

T  F  (  A  (  Y  9  K  )  .  FO  .  0 . 0  )  00  T 0  71 
nn  *  J=KK*M 
A  (  K  9  J)  =A  (K  *  J  )  /A  (  K*K  > 
no  6  T  =  K  K  »  N 
Msfl  (  T  ,<  ) *A ( K  * J ) 

A  (  I  , J ) =  A  (  I  * J ) -W 

y  fr  (  n  A  R  F  (  A  (  f  ,  ,J  )  )  -  .  0  0  0  O  0  1  *0  AP  S  (  W  )  ^  6  1  »  4  <>  4 

4  1  A  (  I  , J ) =  0  . 

^  CONTINUE 
y  =  y  kr 

TF (K.NF.N)  00  TO  1 

70  n=A ( M ,N ) *D 

71  RFTIJRN 
FMO 

FUPPOHTIMF  OTVTDF  (P»MULT»N) 
n o  ( J  p  L  F  PRECISION  d ( 1 i  *11 ) *  M  U  L  T  *  C ( 1  1  1*0(11 ) 
Ml  I  L  T  =  1  • 
no  i  1  =  1  *  N 
r ( T  )  =  p (  I  * i  ) 

T  F ( 4 (  T  )  . FO. o . o  )  0 (  I  )  =1  .  n 
1  OOMT  TNI  IF 
no  o  T  =  1  * n 
no  7  J=  1  *M 

T  P  (P(T,J).FO. 0.4)00  TO  7 
R(T,J)=P.  (T,J)/C(T) 

7  fOMTTM'.lF 

no  o  j= i *n 
n  (  j  )  =  o  (  1  *  J  ) 


IF(n(j)«F0.P,0)D(J)=1 .0 
*  rONTIMUF 
no  4  j=l ?M 
no  a  T  =  ]  ,M 

TF (P ( T , J ) .F0.n,n)00'  TO  4 
D(I,J)=p(T,J)/n(j) 

4  ronjiMiir 

n  F  T!  !P|M 
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